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The availability of ab initio electronic calculations and the concomitant techniques for deriving 
the corresponding lattice dynamics have been profusely used for calculating thermodynamic and 
vibrational properties of semiconductors, as well as their dependence on isotopic masses. The latter 
have been compared with experimental data for elemental and binary semiconductors with different 
isotopic compositions. Here we present theoretical and experimental data for several vibronic and 
thermodynamic properties of CuGaS2 , a canonical ternary semiconductor of the chalcopyrite family. 
Among these properties are the lattice parameters, the phonon dispersion relations and densities 
of states (projected on the Cu, Ga, and S constituents), the specific heat and the volume thermal 
expansion coefficient. The calculations were performed with the ABINIT and VASP codes within the 
LDA approximation for exchange and correlation and the results are compared with data obtained 
on samples with the natural isotope composition for Cu, Ga and S, as well as for isotope enriched 
samples. 

PACS numbers: 63.20.-e, 63.20.dk, 63.20.D-, 68.35.bg, 65.40.Ba, 71.55.Gs, 71.70.Ej 
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I. INTRODUCTION 



The availability of ah initio electronic calculations and 
the concomitant techniques for deriving the correspond- 
ing lattice dynamics have been profusely used in the past 
decade for calculating thermodynamic and vibrational 
properties of semiconductors, as well as their depen- 
dence on isotopic masses. The latter have been compared 
with experimental data for elemental and binary semi- 
conductors with different isotopic compositions Here 
we present theoretical and experimental data for sev- 
eral vibronic and thermodynamic properties of a canon- 
ical ternary semiconductor of the chalcopyrite family: 
CuGaS2. Two main groups of chalcopyrites are usually 
considered, one, denoted as Tin-Vl2, derived from the 
IT VI compounds with zincblende structure, the other, 

II- IV-V2 , derived from the III-V zincblende compounds. 
Examples of the first group are CuGaS2 and AgGaS2, 
whereas the second type is represented, for example, 
by ZnGcAs2. The chalcopyrite structure has the space 
group lA2d and the class 42m with two formula units 
per primitive cell and a longitudinal distortion along the 
c-axis which converts the tetrahedral primitive cell (PC) 
into a tetragonal one. The lattice constants of the tetrag- 
onal PC arc a — b (along x and y) and c (along z)<^ 
The regular tctrahedra with the anion at the center and 
the cations at the vertices are distorted because, e.g., 
of the different lengths of the I-II and the III-II bonds. 
The anion distortion is usually chosen to be along the 
X direction and equal tou-1/4, m = 1/4 corresponds 
to no tetrahedral distortion. In this article we discuss 
lattice parameters and vibrational properties of some T 

III- VI2 chalcopyrites (I = Cu, Ag). The motivation for 
this choice (as opposed to the II-VTV2 materials) is that 



the copper and silver chalcopyrites have received consid- 
erable attention for the production of photovoltaic cells. 
Their energy gaps cover the range 1 - 3.5 eV, i.e., most 
of the frequency of the solar spectrum. From the funda- 
mental point of view, these materials have the property 
that the 3d core electrons of Cu and the id of Ag over- 
lap (and thus hybridize) with the top of their valence 
bands, and thus giving rise to a number of interesting 
anomalies involving negative spin orbit splittings (-0.016 
eV for CuGaS2, (Ref. Q) and nonmonotonic behavior of 
the energy gap versus temperaturci^ 



In this article we focus on the chalcopyrite CuGaS2. 
These sulfides, CuInS2 and the corresponding selenides 
are being considered, together with their alloys, as ef- 
ficient photovoltaic materials. From the fundamental 
point of view CuGaS2 has received less attention than 
the other related chalcopyrites, a reason why here we 
concentrate on the physical properties of this material. 
CuGaS2 was first synthesized by Hahn et al.X They also 
determined by X-ray diffraction the crystal structure and 
the lattice parameters a = h, c, and u of this and 19 other 
related chalcopyrite compounds. Many of this parame- 
ters agree reasonably well with those determined experi- 
mentally up to date and also with recent ab initio calcu- 
lations. CuGaS2 was first found as a mineral (gallite) in 
Namibia and in the Kongo 



CuGaS2 crystallizes with the chalcopyrite structure 
which is closely related to that of zincblende with a slight 
distortion resulting from the tetrahedral bonding of the 
latter ii 
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II. THEORETICAL DETAILS 

The calculations reported here concern the lattice pa- 
rameters of CuGaS2 , the Raman, ir active, and silent k = 
phonons , their Griineisen parameters and the phonon 
dispersion relations, the densities of phonon states (DOS) 
(including the projections on the vibrations of the 3 com- 
ponent atoms), and the optically active densities of two- 
phonon states. In addition, we present ab initio calcula- 
tions of the elastic constants and the bulk moduli Bq and 
Bq' . Because of the large number of phonon at fc = 
(24) we surmised that for a calculation of the volume ex- 
pansion coefficient vs. T a BZ sampling using only the 24 
Griineisen parameters of these phonons (at A: « 0) would 
yield a reasonable approximation to the scarce experi- 
mental results available. This conjecture turned out to 
be correct. Finally, we used the calculated phonon den- 
sity of states (PDOS) to evaluate the specific heat at con- 
stant volume (and the expansion coefficient to evaluate 
the measured constant pressure counterpart) . These cal- 
culations were performed with the natural isotopic abun- 
dance of the constituents of CuGaS2 and also for crystals 
composed of isotopically pure atoms. 

The calculations were based on ab initio electronic 
band structure determinations using density functional 
theory with either the ABINIT or the VASP code3"i^ 
In the ABINIT calculation, normalized pscudo potentials 
were generated by using the Fritz Haber Institute code 
with a valence electron configuration of 3(i^° 4s^ for Cu, 
3dio 4s2 4pi for Ga and 3s^ 3p^ for The wave func- 
tion was expanded in plane wave up to an energy cutoff 
of 40 Ha and the Brillouin zone was sampled by using 
the Monkhorst-Pack method with a 6 x 6 x 6 k-point 
grid. While most electronic calculations were performed 
without spin-orbit (SO) interaction, in order to reveal 
possible effects of this interaction on the lattice proper- 
ties as well as the negative sign of the SO splitting at the 
top of the valence band, a few band structure calculations 
with the VASP code including SO interaction were also 
performed. The effects of this interaction on the lattice 
parameters and dynamics was found to be insignificant. 
We display in Fig. [T]the electronic band structure calcu- 
lated with the VASP code using the LDA exchange- cor- 
relation potential without SO interaction. The details of 
the ab initio electronic band structure calculations have 
been given in Ref. [isl . 

This band structure is similar to those reported re- 
cently by Soni et al. and Brik using a GGA exchange- 
correlation potcntialJ^ii^ Notice that in our case and 
that of Soni et al. the direct gap (1-2 cV) is consid- 
erably smaller than the experimental one (2.43 eV), a 
rather general property of LDA calculations sometimes 
referred to as the " gap problem" Brik brought the cal- 
culated gap to agree with the experimental one by using 
the so-called 'scissors operator'J^ 

An interesting feature of the band structures of Fig. [1] 
and Ref. [l^ is the fact that the bands which correspond 
to the 3d electrons of the copper overlap with the 2p 
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FIG. 1: Electronic band structure of CuGaS2 calculated with 
the VASP code using the LDA exchange-correlation potential 
without SO interaction. The four lowest bands involve mainly 
3s electrons of S. The bands between and -8 eV correspond 
to twenty 3d electrons of Cu and twenty four 3p of S. Notice 
that the calculated energy gap (1.1 eV) is much smaller than 
the experimental one (2.4 eV, Ref. 0), reflecting the so-called 
"gap-problem".— The notation of the special points in the 
Brillouin zone (T - T - N - P) is identical to that of Ref. [l3. 



bands of sulfur. This peculiarity results in a negative SO 
splitting at the top of the valence bands (-0.016 eV) (ReL 
|20[ ). an anomalous sign also observed in CuCl (ReL [2l[ ). 
ZnO (ReL HI), and /3-HgS (ReL Here we shall 

no longer discuss electronic properties of chalcopyrites 
and concentrate on the lattice parameters, phonons and 
thermodynamic properties. 

The reader may wonder why we use two different DFT 
codes. The reason is that we are familiar with both codes 
and, while we know that they lead to similar results in the 
case of monoatomic and binary crystals, it is not obvi- 
ous that they also will do so for more complicated struc- 
tures. This was shown to be so in the case of cinnabar 
(a- HgS^), with three atoms per PC. Here we examine 
chalcopyrites, with 8 atoms per PC. 

The wave function representation for the two codes is 
different as well as the methodology e. g. to calculate the 
phonon vibrational spectra. Therefore, the use of both 
codes is complementary. In some cases, the ABINIT code 
is more computationally demanding but the precision can 
be increased. Therefore, for some properties, such as the 
phonon spectra, we have employed the ABINIT imple- 
mentation where density functional perturbation theory 
was usedi ^^i^^ 



III. EXPERIMENTAL DETAILS 

The CuGaS2 crystals investigated here were grown by 
a vapor phase transport technique using iodine as trans- 
port agent. The isotopically nearly pure (99.5% ^"'S, 
99.9% ^^Cu and 99.6% ''^Ga) elements were purchased 
from Trace Science International, Ontario, Canada. 
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For the single crystal X-ray structure determination 
CuGaS2 crystals were selected under the polarization mi- 
croscope for crystallographic investigations. As many of 
the crystals showed twinning according to 1, with 1 be- 
ing the identity operation, a small cuboid fraction (0.24 
X 0.20 X 0.14 mm'^) of a yellow transparent crystal was 
oriented on a four-circle diffractomcter using graphite- 
monochromatizcd Ag-Ka radiation (CAD 4, Enraf No- 
nius, Delft, The Netherlands). The lattice parameters 
of the tetragonal lattice (space group I42d, No. 122) 
were refined from 25 centered high-indexed reflections to 
a = 5.3512(6) and c = 10.478(3) A. Bragg intensities in 
three octants of the Ewald sphere (-9 < ii, -k < 9, -18 
< J < 18, 3.4° < i9 < 29.9°, R„t = 0.0648 after merg- 
ing) were collected in the cj-2z9 data acquisition mode and 
corrected for Lorentz, polarization and absorption effects. 
The structure refinement was performed with full matrix 
least-squares cycles on {Rl = 0.0472 for I > 2 al and 
0.0832 for ah I, wR2 = 0.1178 for J > 2 crJ and 0.1350 for 
all)^ 

Temperature dependent lattice parameters and the 
thermal expansion coefficients were determined by pow- 
der X-ray diffraction (A = 1.54 A) on crushed crystals 
grown from elements with the natural isotopic abun- 
dance. 

The heat capacities were measured on samples of typi- 
cally ~ 20 mg between 2 and 280 K with a physical prop- 
erty measurement system (Quantum Design, San Diego, 
CA) as described in detail in Ref. [2^. Between room 
temperature and 1100 K the heat capacities of a ~100 
mg polycrystalline sample were determined with a DSC 
404 Fl Pegasus differential scanning calorimeter (heat- 
ing rate 20 K/min) with the sample kept in an argon 
atmospherei^ Up to 1100 K a reduction of the sample 
mass was not observed. 



IV. RESULTS AND DISCUSSION 
A. Crystal Structure 

The available structural parameters of CuGaS2 exhibit 
considerable dispersion, especially the x positional pa- 
rameter of the S atoms which is very close to 1 /4 (in the 
forthcoming called u) has so far been determined with 
limited reliability and found to be larger than 1/4. Ad- 
ditionally, u has been found to vary non-monotonically 
from CuAlS2, via CuGaS2 to CuInS2i^ In order to in- 
crease the accuracy of u we redetermined the crystal 
structure of CuGaS2 using high-quality single crystals 
and up-to-date X-ray diffraction techniques which allow 
us to decrease the experimental error in w by a factor 
of ~7. Table ITV^ summarizes the results of our crystal 
structure redetermination. 

The positional parameter u is clearly smaller than 1 /4 
-2.5%) and decreases monotonically throughout the 
series CuAlS2 - CuGaS2 - CuInS2. Table HIl summarizes 
the experimental and calculated structural parameters of 



CuXS2 (X = Al, Ga, In) as obtained from our calcula- 
tions and as available in the literature. 

Since the available calculations did not include SO 
splitting we performed, for comparison, calculations with 
SO splitting. They are displayed in Table IIIII together 
with a VASP-NOSO calculation carried out by Chen 
et al^ using a GGA DFT for the exchange-correlation 
hamiltonian. 



B. Elastic Properties 

The ABINIT code used for our calculations provides 
also the stiffness constants, Cijr^ six independent ones 
in the case of the chalcopyrites (Table ITV)) . 

We have not been able to retrieve experimental values 
for these from the literature. In Table IIVI we display 
for comparison data for CuGaSe2 and AgGaS2. For the 
latter experimental data arc also quoted. We have also 
included bulk modulii, Bq, and their pressure derivative, 
Bq', for all three compounds. Except for Bq', all elastic 
parameters undergo a monotonic decrease through the 
series CuGaS2, CuGaSe2 to AgGaS2, which is likely to 
be related to the corresponding increase of the lattice 
parameters. 
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FIG. 2: (color online) Phonon dispersion relations of CuGaS2 
as calculated with the ABINIT-LDA code. The notation of 
the special points in the Brillouin zone (T - F - N - P) is 
identical to that of Ref. [l^. The (red) dots represent some of 
the available experimental frequencies (see Table . 

The calculated phonon dispersion relations are shown 
in Fig. [2] in the reduced Brillouin zone. The F - T direc- 
tion corresponds to [001] whereas F - N corresponds to 
[110] For comparison, wc have added a few points ob- 
tained from Raman and infrared spectroscopy measure- 
ments given in more detail in Table |V] where they are 
also compared with our ab initio calculations and those 
of Akdogan et al.M-^ 

Figure [3] displays the phonon densities of states corre- 
sponding to the motion of the three constituent atoms 
calculated from the dispersion relations shown in Fig. [21 
As expected, the low- frequency band - 100 cm~^ cor- 
responds mainly to Cu and Ga vibrations whereas the S- 
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TABLE I: Summarized results of the crystal structure determination (T = 293(2) K) of CuGaS2 (space group /42rf, No. 122). 
Standard deviations of the last digit are given in parentheses. The lattice parameters were determined to be a = 5.3512(6) A, 
c — 10.478(3) A resulting in a unit cell volume of 300.03(9) A'^ which contains 4 formula units {Z= 4). Further experimental 
details are deposited under the no. CSD-422615 at the Fachinformationszentrum Karlsruhe.™ 

Fractional atomic coordinates and equivalent isotropic displacement parameters (A'') 
Atom Wyckoff-Nr. x y z (7cquiv 

Ga 4a 0.0078(4) 

Cu 4b 5 0.0159(4) 

S M 0.2437(5) \ | 0.0085(4) 

Anisotropic displacement parameters Uij{K'') 
Atom Uii U22 U33 U12 U13 U23 

Ga 0.0080(4) 0.0080(4) 0.0076(5) 
Cu 0.0155(5) 0.0155(5) 0.0165(6) 

S 0.0113(11) 0.0051(10) 0.0092(6) -0.0016(5) 



like contributions are mainly above the ^280 cm~^-gap. 
we note that there are also some S-like contributions be- 
low '^120 cm~^ originating from Cu - S vibrations. The 
partial densities of states are e.g. useful for calculating 
the effect of isotope disorder on the phonon linewidths.— 
We have not found in the literature second-order Ra- 
man spectra of CuGaS2 which would correspond to the 
sum and difference spectra of FigH) Nevertheless, it is 
possible to establish a correspondence between the calcu- 
lated two-phonon Raman spectra of CuGaS2 and (Fig|3]) 
and the measured ones of /3-ZnS shown in Fig. 1 of 
Ref. HO- We present here the calculated sum and dif- 
ference densities of states of CuGaS2 in the hope that 
they will help in interpreting measured spectra when they 
become available for the CuGaS2 or other chalcopyrite 
compounds. 

Table |V] also lists the calculated and experimental 
Griineisen parameters of F-point phonons. As usual, 
most of the the Griineisen parameters are positive ex- 
cept some at the lowest frequencies. 

The effect of lattice vibrations on the volume Vb(T) of 
a (cubic) crystal can be expressed in terms of the mode 
Griineisen parameters 7qj and the mode frequency as 



AFo(r) 

Vo 



BoV 



^13 ^13 ^^10) + |] ' (1) 



where is the Bose-Einstein factor 



TABLE III: Lattice parameters, fractional x atomic coordi- 
nate of sulphur (u), cell volume and bulk moduli, Bo and 
Bq as obtained from our ah initio calculations either with 
spin-orbit coupling (SO) or without (NOSO). 



code 


a (A) 


c(A) 


u 


Vcoii (A") 


Bo (GPa) 


Bo' 


ABINIT-NOSO 


5.262 


10.45197 


0.244 


289.4 


93.3 


4.7 


VASP-NOSO-LDA 


5.226 


10.380 


0.2450 


283.5 


94.3 


4.5 


VASP-SO-LDA 


5.2258 


10.3818 


0.2450 


283.5 


91.9 


5.1 


VASP-NOSO-GGA° 


5.3700 


10.643 


0.2491 


306.9 


85 


4.7 



Ref. 30 



TABLE IV: Comparison of the stiffness constants obtained 
from ABINIT calculations performed for the three isostruc- 
tural compounds CuGaS2, CuGaSe2 and AgGaS2 with exper- 
imental data obtained in this work and in the literature. 

parameter CuGaS2 
calc exp 



Cii (GPa) 
C12 (GPa) 
Ci3 (GPa) 
C33 (GPa) 
C44 (GPa) 
C66 (GPa) 
Bo 

Bo' (GPa) 



132.23 
78.4 
79.4 
144.1 
56.1 
56.3 
96.4 
4.5 



94^= 



97" 



6.3^*, 4^= 



CuGaSe2 


AgGaS2 


calc 


exp 


calc 


exp 


112.2° 




85.3" 


86.5'' 


66.4° 




52.4'' 


56.0* 


68.1° 




59.9'' 


59.6* 


113.2° 




76.2'' 


75.1* 


48.4° 




32.4'' 


24.9* 


48.5° 




36.1'' 


31.4* 


83° 


102° 60"*, 62.3'' 


72.2* 






4.2'' 


4* 



° Ref. 32 
* Ref. 33: 
" Ref. 34 
Ref. 35 
" Ref. 36 



average values of references therein. 



TABLE II: Comparison of the results of VASP-NOSO-LDA 
calculations performed for the three isostructural compounds 
CuXS2 (X = Al, Ga, In) with experimental data obtained in 
this work and in the literature. 



parameter CUAIS2 



CuGaS2 



CuInSa 



VASP exp VASP exp VASP exp 

a = 6 (A) 5.2055 5.326 5.226 5.3512(6) 5.482 5.5221 

c (A) 10.3765 10.436 10.380 10.478(3) 10.9301 11.1043 

u 0.2505 0.271 0.2450 0.2437(5) 0.21651 0.2145 

° Ref. 2^ and average values of references therein. 



(2) 



and V and Bq are the volume and the bulk modulus 
of the crystal, respectively.— 

Because of the large number of phonon bands we sur- 
mised that the Griineisen parameters of the phonons 



5 




100 200 300 

E (cm"') 



400 



FIG. 3: (color online) Phonon density of states (PDOS) pro- 
jected on the three constituent atoms. Note the gap between 
the contributions of (essentially) Cu and Ga (below 200 cm~^) 
and those above ~ 280 cm^'^ which are essentially S like. 
There is considerable S like weight below ~120 cm~^. The 
sum of these partial densities of states has been used to cal- 
culate the temperature dependence of heat capacities. 
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FIG. 4: Sum (upper panel) and difference (lower panel) 
phonon densities of states of CuGaS2 as calculated from the 
dispersion relations shown in Fig. |4l 



TABLE V: Calculated and measured phonon frequencies of 
CuGaS2 (in cm~^) at the center of the BZ. Also, averaged 
experimental values reported in about 9 different publications 
(see Akdogan et al^, Table 4). The values of the correspond- 
ing Griineisen parameters, 7, calculated by us are also listed 
and compared with measurements. 
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at F suffice, to a first approximation, for tfie evalua- 
tion of tlic volume thermal expansion coefficient, av= 
(l/Vg) dVo{T)/dT, to include in the summation in Eq. 
([5]) the phonons and their Griineisen parameters at the 
F-point, where Vq represents the volume of the primi- 
tive cell. As usual in semiconductors, negative values 
of "fqj lead to negative thermal expansion coefficients at 
low temperatures4^ The temperature dependence of Vq 
is given by 



AVojT) 
Vo 



BnV 



Wj[ns(wj) -f 



(3) 



where uj and jj arc the zone-center phonon frequencies 
and the Griineisen parameter, respectively. 

The volume thermal expansion coefficient can also be 
obtained from the variation of the entropy, S{P, T), with 
pressure via the thermodynamic relationship 



1 f dS{P,T) 



(4) 
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We calculated within the ABINIT code the entropy at 
ambient pressure and for pressures of 0.4 and 1 GPa and 
took the numerical derivatives. In Fig. [5] wc compile 
these results and the the results of our calculations using 
the F-point phonon frequencies and the mode Griineisen 
parameters summarized in Table |V] (and its temperature 
derivative), i.e. the volume thermal expansion coeffi- 
cient, ay(r) = (1/Vb) dV{T)/dT, with literature data 
by Bodnar et ai, Schorr et ai, and with our high tem- 
perature X-ray data. For comparison we also display the 
volume thermal expansion coefficient of Reeber and Pow- 
eh et al.^d^ 



pacity experiments (sec below). The temperature de- 
pendence of the volume thermal expansion coefficient for 
CuGaS2 is expected to be similar to that of the iso- 
baric and isoelectronic ZnS. This is indeed the case; up 
to room temperature our calculations coincide remark- 
able well with volume thermal expansion data for ZnS 
obtained by Reeber and Powell.— Above room temper- 
ature Roberts et al. observed a linear increase of the 
volume thermal expansion coefficient of zincblende ZnS 
at a rate of 1.16xlO~^K~^ which is not found in our cal- 
culations for CuGaS2'^ This linear increase of the vol- 
ume thermal expansion coefficient of wurtzite ZnS and 
zincblende ZnS has recently also been obtained by ab 
initio calculations.— 
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C. Heat Capacity 

The phonon density of states displayed in Fig. [3] allows 
calculation of the free energy F{T) by using the expres- 
sion: 



F{T) 



(■^ + kBTln[2nB{io)])p{uj)(hj (5) 
^ 



and of the specific heat at constant volume by taking 
the second derivative of the free energy 



FIG. 5: (color online) (red) dotted line: Temperature depen- 
dence of the volume thermal expansion coefficient of CuGaS2 
as obtained from the zone-center Griineisen parameters 7 and 
frequencies as given in TableFVl using eq. ([3}. (black) solid line 
thermal expansion coefficient from the pressure dependence of 
the entropy (eq. (0])). Literature data by Reeber et al. and 
Roberts et al. and our X-ray results are also given.— 

The volume thermal expansion coefficients obtained 
from the F-point phonon frequencies and from the 
pressure-derivative of the entropy are in good agree- 
ment, thus justifying the approximation used in eq. ([3]). 
The results of our calculations clearly reveal a nega- 
tive thermal volume expansion coefficient at low tem- 
peratures as has been found experimentally by Schorr 
and Scheptyakov4i However, at high temperatures our 
calculated data deviate markedly from the experimen- 
tal findings, they are about 15% lower than the data 
by Bodnar and Orlova obtained by X-ray diffraction 
experiments!^ However, there still appears to be some 
scatter in the experimental data. While Bodnar et al. 
report a room temperature volume thermal expansion 
coefficient of 2AxlQ-^K-^ close to 26.5x10-6K-i ob- 
tained by Yamamoto et al., Malsagov et al. found a 
room temperature value of only 19.2x 10~^K~^, close to 
the results of our calculations Our calculations 
also agree rather well with our high-temperature X-ray 
diffraction data carried out on a polycrystalline sample 
prepared from the same crystals used for the heat ca- 



In Eq. ([S]), fcs is the Boltzmann constant, ns the Bose- 
Einstein factor, and p{uj) the phonon density of states. 
The high frequency cutoff of the latter defines the upper 
limit of integration in Eq. ([S]). 

The difference between the calculated heat capacity at 
constant volume, Cy, and that at constant pressure, Cp, 
(the quantity that is experimentally obtained) is related 
to the volume thermal expansion coefficient, ay, accord- 
ing to 

Cp (T) - C'v (T) ^al{T)-B- V^oi ■ T, (7) 

where Bq is the (isothermal) bulk modulus and Vmoi 
the molar volume. With av{T > 300K) « 2.3 x 
10~^K~^ the contribution to the heat capacity from ther- 
mal expansion amounts to 0.7 J/molK and 2.3 J/molK 
at 300 K and 1000 K, respectively. 

In the literature heat capacity data for CuGaS2 are 
available in the low temperature regime (13 K < T < 38 
K) and above room temperature up to ~600 Ki^— Our 
data connect these temperature regimes and extend the 
temperature range up to 1100 K. At low temperature our 
results are in good agreement with the data of Abrahams 
and Hsu with improved resolution of the maximum in 
Cp/T^ (cf. Fig. [T]), at high temperatures they connect 
well to the data by Neumann et ai. In Fig. [5] we display 
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the heat capacity of CuGaS2 over a temperature range 
from 2 K up to 1100 K together with the theoretical 
results based on the PDOS shown in Fig. [3l To account 
for the contributions of the thermal expansion, becoming 
especially noticeable above room temperature, we have 
used our data of the volume thermal expansion coefHcient 
extended by a linear increase above room temperature of 
1.16xlO~^K~^. identical to that found by Roberts et al. 
for zincblendc ZnS4^ as displayed in Fig. [5] and used Eq. 
([7]) to calculate Cp. 
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FIG. 6: (color online) Temperature dependence of the mo- 
lar heat capacity of CuGaS2. (o) and (blue) dashed line: our 
experimental data, (red) dotted line: results of ABINIT calcu- 
lations of Cv, (black) solid line: Cp according to Eq. (O us- 
ing our calculated volume thermal expansion coefficient oi{T) 
displayed in Fig. (5] Above room temperature we have added 
a linear increase of alT) which amounted to 1.16xlO~^K~^, 
identical to that found by Roberts et al. for zincblende ZnS.— 
Literature data by Abrahams and Hsu and Neumann et al. 
are also displayed.— The vertical line indicates the Petit- 
Dulong value of 12x7?, where r is the molar gas constant. 

In some of our previous worksi^iS we have investigated 
the dependence of Cp/T'^ on the isotopic masses of the 
constituents of the compounds (up to binary compounds, 
so far) and compared the experimental data with the the- 
oretical results. For elementary and binary compounds 
we have also investigated the relationship of the logarith- 
mic derivatives versus temperature and versus the masses 
of the constituentsi^ 

In all logarithmic derivatives we observe a peak cen- 
tered at about 20 K. For S there is an additional broad 
band with its maximum at about '^100 K. The maxima 
in the logarithmic derivatives reflect the structure of the 
PDOS projected on the corresponding atoms. Clearly, 
the broad high energy feature visible exclusively for S 
originates from the PDOS with mainly S character be- 
tween 280 and 380 cm~^. The ratio of the maximum 



100 



FIG. 7: (color online) Temperature dependence of the molar 
heat capacity of CuGaS2. (o): our experimental data. Liter- 
ature data by Abrahams and Hsu are also given.— 



frequency (converted into Kelvin) in the PDOS and the 
maximum temperature of the logarithmic derivative is 
~6, a ratio which has been found in a number of previ- 
ous investigations^^ Cu shows the highest feature in the 
logarithmic derivative at low temperature reflecting the 
2g(^igh weight of the Cu projected PDOS at low energies. 
We have demonstrated that there is a close relation- 
ship of the logarithmic derivatives of the heat capaci- 
ties with respect to temperature and isotope mass<^ The 
straightforward extension of the relationship the logarith- 
mic derivatives versus temperature and versus the masses 
of the isotopes to the ternary compound, CuGaS2 in our 
case, is given by 



1 rfin(cyr3) 

2 ^-^^ d\nT > 



d\n{Cp/T^) d\n{Cp/T^) d\n{Cp/T^) 



dlnMcu 



dlnMca 



din Ms 
(8) 

where Mqu, Afca, and Ms arc the masses of the 
three constituents, i.e. Cu, Ga, and S, respectively, for 
CuGaSa. 

Figure [5] conflrms the relationship of the temperature 
dependence of the logarithmic derivatives with respect to 
temperature and to the isotope masses established ear- 
lier by us is also valid for multinary (in the present case 
ternary) compounds. 

We have demonstrated that for low temperatures, T — ^ 
0, the logarithmic derivatives are related to the ratios of 
the atomic mass to the molar mass according to^ 



dlnCjT^ _ 3 Mcu 
rflnMcu ~ 2 Mcu + Mca 

dlnq/T^ ^ 3 A/Ga 
dlnMca 2 Mcu + Mca 

dlnC„/T3 _ 3 2Ms 
d In Ms ^ 



2Ms 

= 0.53 



2Ms 



2 Mcu + A'/Ga + 2A'/s 



= 0.48 (9) 
(10) 
(11) 



0.49. 
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In the case of CuGaS2 these three ratios arc, fortu- 
itously, approximately equal to 0.5. 



V. CONCLUSIONS 



Ab initio electronic band-structure techniques, espe- 
cially those which use up-to-date computer codes like 
VASP or ABINIT, are powerful methods to investigate 
electronic, optical, vibronic and thermodynamic proper- 
ties of crystals and the results with experimental data. 
Here we apply these techniques to CuGaS2 which has 
chalcopyrite structure (space group 142(1 (No. 122), two 
molecules per primitive cell.) more complicated than 
those usually dealt with. We used the ABINIT code 
to calculate the frequencies of Raman and IR phonons 
and their dispersion relations. The densities of states of 
one and two phonons have also been calculated. We de- 
vote the last section to present experimental data on the 
specific heat versus temperature of samples grown with 
the natural isotopic abundances and those grown with 
isotopically modified ones. These results are compared 
with ab initio calculations. Generally, good agreement 
between experiment and ab initio results is obtained. 
A redetermination of the crystal structure parameters 
is presented which decreases the discrepancies with ex- 
perimental data apparent in the literature. The peaks of 
the Cu and Ga mass derivatives in Fig. |8] are at simi- 
lar temperatures. In the process of measuring AgGaS2 
wherein the masses of the two cations are considerably 
different we observe a clear difference of the low temper- 
ature phonon spectrum related to the differences of the 
atomic masses of Ag and Ga. 



100 
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FIG. 8: (color online) (a - c) Logarithmic derivatives of the 
specific heats, Cp/T^, with respect to the atomic masses of the 
constituents, Cu, Ga and S (from top to bottom), respectively. 
The (black) solid lines represent the logarithmic derivatives 
calculated from our theoretical heat capacities (ABINIT). (d) 
(red) o Logarithmic derivatives of the specific heat of CuGaS2 
(natural isotope composition) with respect to the tempera- 
ture. Plotted is the quantity i x {dlnCp/T'-^ /dlnT+3) which 
is compared with the sum of the logarithmic derivatives with 
respect to the isotopic masses (blue) □: experiment; black 
solid line: sum of theory curves shown in (a - c). At very 
low temperatures the theoretical logarithmic derivatives have 
been extrapolated (dashed lines) to the values imposed for 
T ^ by eq. dUl}. 



Acknowledgments 

AHR has been supported by CONACYT Mex- 
ico under projects J-83247-F, Binational Collabora- 
tion FNRS-Bclgium-CONACYT and PPPROALMEX- 
DAAD-CONACYT. A. M. acknowledges the finan- 
cial support from the Spanish MCYT under grants 
MAT2010-21270-C04-03, CSD2007-00045 and the super- 
computer resources provides by the Red Espahola de 
Supercomputacion. We are grateful to Y. Pouillon and 
O. Castillo for valuable technical and computational sup- 
port. We are also indebted to E. Kotomin for a critical 
reading of the manuscript. 



9 



* Corresponding author: E-mail [R.Kremer@ fkf.mpg.de' 

^ M. Sanati, S. K. Estreiclier, and M. Cardona, Solid State 

Commun. 131, 229 (2004) 
^ A. Gibin, G.G. Devyatykh, A.V. Gusev, R.K. Kremer, 

M. Cardona, and H.-J. Pohl, Sol. State Comm. 133 569 

(2005). 

^ M. Cardona, R. K. Kremer, R. Lauck, G. Siegle, A. Munoz, 
and A. H. Romero, and A. Schindler, Phys. Rev. B 81 
075207 (2010). 

* J.E. Jaffe and A. Zunger, Phys. Rev. B 28, 5822 (1983). 

^ H. Horinaka, Y, Yamamoto, and T. Miyauchi, Jap. J. Appl. 

Phys 17, 521 (1978). 
® J.Serrano, Ch. Schweitzer, C.T. Lin, K. Reimann, 

M.Cardona, and D. Frohlich, Phys. Rev. B 65, 125110 

(2002). 

H. Hahn, G. Frank, W. Klingier, A. D. Meyer, G. Stoerger, 

Z. Anorg. AUg. Chem. 271, 153 (1953). 
® H. Strunz, B. H. Geier, E. and Seeliger, Neues Jahrbuch 

fiir Mineralogie, Monatshefte 241 (1958). 
^ H. Strunz, B. H. Geier and E. Seliger, Am. Mineral. 44, 

906 (1959). 

1° X. Gonze et al, Comput. Mater. Sci. 25, 478 (2002). 
ABINIT is a collaborative project of the Universite 
Catholique de Louvain, Corning Incorporated and other 
contributors. 

" X. Gonze, B. Amadon, P.-M. Anglade, J.-M. Beuken, F. 
Bottin, P. Boulanger, F. Bruneval, D. Caliste, R. Cara- 
cas, M. Cote, T. Deutsch, L. Genovese, Ph. Ghosez, M. 
Giantomassi, S. Goedecker, D. R. Hamann, P. Hermet, F. 
JoUet, G. Jomard, S. Leroux, M. Mancini, S. Mazevet, M. 
J. T. Oliveira, G. Onida, Y. Pouillon, T. Rangel, G.-M. 
Rignanese, D. Sangall, R. Shaltaf, M. Torrent, M. J. Ver- 
straete, G. Zerah, and J. W. Zwanziger, Computer Phys. 
Commun. 180, 2582 (2009). 

G.Kresse and J. Furthmiiller, Comp. Mat Sci 6, 15 (1996). 
G.Kresse, D. K. Joubert, Phys Rev. B 59, 1758 (1999) and 
references therein. 

M. Fuchs, M. Scheffler, Comput. Phys. Commun. 119, 67 
(1999). 

M. Cardona, R. K. Kremer, G. Siegle, A. Munoz, A. H. 
Romero, and M. Schmidt, Phys. Rev B 82, 085210 (2010). 
^® R. W. Godby, M. Schliiter, and L. J. Sham , Phys. Rev B 
37, 10159 (1988). 

A. Shileika, Surf. Science 37, 730 (1973). 

A. Soni, V. Gupta, C. M. Arora, A. Dashora and B. L. 

Ahuja, Solar Energy 84, 1481 (2010). 

M. G. Brik, J. Phys.: Condens. Matter 21, 485502 (2009). 
^° B.TeU and P.M. Bridengaugh, Phys Rev. B 12, 3330 
(1975). 

K. Shindo, A. Morita, and H. Kamimura, J. Phys. Soc. 
Japan 20, 2054 (1945). 
^2 P. Carrier and Su-Huai Wei, Phys. Rev. B 70, 035212 
(2004) 

M. Cardona, R. K. Kremer, R. Lauck, G. Siegle, A. Munoz, 
and A. H. Romero, Phys. Rev. B 80, 195204 (2009). 
X. Gonze, Phys. Rev. B 55, 10337 (1997). 
X. Gonze and C. Lee, Phys. Rev. B 55, 10355 (1997). 
^® Further information about data collection, data reduc- 
tion, refinement and geometrical details can be ob- 
tained from Fachinformationszentrum Karlsruhe, D-76344 



Eggenstein-Leopoldshafen (Fax: -1-49-7247-808-666, e- 
mail: crysdata@fiz-karlsruhe.de) by quoting the deposition 
number CSD-422615. 

NETZSCH-Geratebau GmbH, Wittelsbacherstr. 42, D- 
95100 Selb, Germany. 

J. Serrano, R. K. Kremer, M. Cardona, G. Siegle, A. H. 
Romero, and R. Lauck, Phys. Rev. B 73, 094303 (2006). 
J. E Jaffe and A. Zunger, Phys. Rev. 28, 5822 (1983). 
3° Shiyou Chen, X. G. Gong, and Su-Huai Wei, Phys. Rev B 
75, 205209 (2007). 

D. R. Hamann, X. Wu, K. M. Rabe, and D. Vanderbilt, 
Phys. Rev. B 71, 035117 (2005). 

C. Parlak and R. Eryigit, Phys. Rev. B 73, 245217 (2006). 
Hon Hai-Jun, Zhu Shi-Fu, Zhao Bei-Jun, Yu You, and Xie 
Lin-Hua, Phys. Scr. 82, 055601 (2010). 
M. Bettini and W. B. Holzapfel, Sohd State Commun., 16, 
27 (1975). 

A. Werner, H. D. Hochheimer, and A. Jayaraman, Phys. 
Rev. B23, 3836 (1981). 

T. Tinoco, J.P. Itie, A. Pohan, A. San Miguel, E. Moya, 
P. Grima, J. Gonzalez, and F. Gonzalez, J. Physique IV C 
9, 151 (1994). 

C. Carlone, D. Olego, A. Jayaraman, and M. Cardona, 
Phys. Rev. 22, 3877 (1980). 

M. Akdogan and R. R. Eryigit, J. Phys.: Condens. Matter 

14, 7493 (2002). 
^® J. Serrano, A. H. Romero, F. J. Manjon, R. Lauck, M. 

Cardona, and A. Rubio, Phys. Rev B 69, 094306 (2004). 

J. Serrano, A. Cantarero, M. Cardona, N. Garro, R. Lauck, 

R. E. Tallman, T. M. Ritter, and B. A. Weinstein, Phys. 

Rev. B 69, 014301 (2004). 
'^^ M. Cardona and M. L. W. Thewald, Rev. Mod. Phys. 77, 

1173 (2005). 

A. Debernardi and M. Cardona, Phys. Rev. B 54, 11305 
(1996). 

I. V. Bodnar and N. S. Orlova, phys. stat. sol. (a) 78, K59 
(1983). 

S. Schorr and D. Sheptyakov, J. Phys.: Cond. Matter 20 
104245 (2008). 

N. Yamamoto, H. Horinaka, and T. Miyauchi, Jap. J. Appl. 
Phys. 18, 255 (1979). 
""^ A. U. Malsagov, I. M. Berfirer, and B. S. Kulbezhev, Izv. 
Vyssh. Ucheb. Zav., Fizika 29, 126 (1986); cited from H. 
Neumann, Cryst. Res. Technol. 22, 723 (1987). 
R. R. Reeber and G. W. Powell, J. Appl. Phys., 38, 1531 
(1967); R. R. Reeber, phys. stat. solidi a 32, 321 (1975). 
R. B. Roberts, G. K. White, and T. M. Sabine, Aust. J. 
Phys. 34, 701 (1981). 

S. Q. Wang, Appl. Phys. Lett. 88, 061902 (2006). 
'^^ S. C. Abrahams and F. S. L.Hsu J. Chem. Phys. 63, 1162 
(1975). 

H. Neumann, G. Kiihn, and W. MoUer, phys. stat. solidi 
B 144, 565 (1987). 

A. H. Romero, M. Cardona, R. K. Kremer, R. Lauck, G. 
Siegle, J. Serrano, and X. C. Gonze, Phys. Rev. B 78, 
224302 (2008). 

M. Cardona, R. K. Kremer, R. Lauck, G. Siegle, J. Serrano, 
and A. H. Romero, Phys. Rev. B 76, 075211 (2007). 



